Signature of odd-frequency pairing correlations induced by a magnetic interface 
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We investigate the mutual proximity effect in a normal metal contacted to a superconductor through a magnetic 
interface. Analytical and self-consistent numerical results are presented, and we consider both the diffusive and 
ballistic regimes. We focus on the density of states in both the normal and superconducting region, and find 
that the presence of spin-dependent phase-shifts occurring at the interface qualitatively modifies the density of 
states. In particular, we find that the proximity-induced pairing amplitudes in the normal metal region undergo a 
conversion at the Fermi level from pure even-frequency to odd-frequency. Above a critical value of the interface 
spin-polarization (or, equivalently, for fixed interface spin-polarization, above a critical interface resistance), 
only odd frequency correlations remain. This is accompanied by the replacement of the familiar proximity 
minigap or pseudogap in the normal layer by an enhancement of the density of states above its normal state 
value for energies near the chemical potential. The robustness of this effect towards inelastic scattering, impurity 
scattering, and the depletion of the superconducting order parameter close to the interface is investigated. We 
also study the inverse proximity effect in the diffusive limit. We find that the above-mentioned conversion 
persists also for thin superconducting layers comparable in size to the superconducting coherence length ^s, as 
long as the inverse proximity effect is relatively weak. Concomitantly, we find a shift in the critical interface 
resistance where the pairing conversion occurs. Our findings suggest a robust and simple method for producing 
purely odd-frequency superconducting correlations, that can be tested experimentally. 



PACS numbers: 

I. INTRODUCTION 

The proximity effect in hybrid structures with superconduc- 
tors offers an arena of interesting physics to explore which 
could also prove to be useful in nanotechnological devices. 
The incorporation of ferromagnetic elements in such hybiid 
structures activates the spin degree of freedom, which has 
a number of important consequences for how the proximity 
effect is manifested in physical quantities.EH in the case of 
a ferromagnet I superconductor (F|S) bilayer, it is known that 
so-called odd-frequency pairing is generated.'* Odd-frequency 
pairing has been studied previously ^ ^ in particular in connec- 
tion with the search for exotic superconducting states that may 
arise via the mechanism of spontaneous symmetry breaking. 
A particular feature of such odd-frequency paring states is a 
strong retardation effect which makes the equal-time coiTela- 
tor vanish for the Cooper pair 

Apart from the mechanism of spontaneous symmetry 
breaking, odd-frequency pairing correlations can also be cre- 
ated by an induced symmetry breaking. The general re- 
quirement for such a generation of odd-frequency correla- 
tions is that either translational symmetry (for odd-frequency 
singIet),'2E^E^ or both translation al an d spin-rotational sym- 
metry (for odd-frequency triplet), are explicitely bro- 
ken. As a result, one would expect to see odd-frequency su- 
perconductivity as a quite generic feature of proximity struc- 
tures. Although this fact is well known since long among the 



community dealing with inhomogeneous problems in super- 
conductivity, it is only recently that the attention has shifted 
to the question: how may one extract and detect these ex- 
otic pairing correlations, and in particular the odd-frequency 
triplet state, experimentally? 

There are two major difficulties associated with the de- 
tection of the odd-frequency triplet state. One obstacle is 
that such a state induced in F|S bilayers often has a very 
short penetration depth into the ferromagnetic region of or- 
der ^ Oinm). In fact, unless there are magnetic inhomo- 
geneities present in the interface region,G21 it is limited by the 
magnetic coherence length which usually is much smaller 
than the superconducting coherence length ^s- A second ob- 
stacle related to the detection of odd-frequency correlations is 
that these often compete with even-frequency superconduct- 
ing correlations in the same mateiial, masking their presence. 
To find smoking gun signatures of odd-frequency pairing is 
therefore a rather challenging issue to tackle, although there 
are a few experimental works whic h hav e pointed towards fin- 
gerprints of odd-frequency pairing.CSIll 

Recently, it has been realized that the interface properties in 
hybrid structures with superconductors play a pivotal role in 
magnetic aspects of the proximity effect.'^^EIl In most works, 
non-magnetic (or spin-inactive) interfaces have been consid- 
ered, even in the presence of feiTomagnetic elements. Utiliz- 
ing the quasiclassical theory of superconductivity, such inter- 
faces are modeled as effective boundary conditions. For the 
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FIG. 1; (Color online) Proposed experimental setup for ob- 
servation of the odd-frequency component in a normal metal 
layer| superconductor junction. 



general Eilenberger equation, boundary conditions for non- 
magnetic systems and spin-inactive interfaces were first for- 
mulated in implicit form in Refs .1221 and l23l An explicit for- 
mulation has been derived in Ref.|24l For the diffusive limit of 
the theory, described by the Usadel equation, boundary con- 
ditions have been formulated in Refs.l25]andl26l 

However, the spin-dependent properties of the interface 
may become important when ferromagnetic elements are 
present in the system. In particular, the transmission proper- 
ties of spin-t and spin-^ electrons into a ferromagnetic metal 
are different, which gives rise to both spin-dependent con- 
ductivities (spin filteringj^S a nd spin -dependent phase shifts 
(spin-DIPS) at the interface .tHEliKlS! a generalization of 
boundary conditions to spin-active interfaces was given in 
Refs. [35l - B7l that has been generalized to include systems with 
strong exchange splitting of the energy bands in Refs. TT' and 

m 

The spin-DIPS can lead to qualitatively novel effects in su- 
perconducting hybrid systems. Very recently, the proximity 
effect in a normal metal [superconductor (N|S) bilayer with a 
magnetic interface was studied in Ref. 39, and a surprising 
result was unveiled. Namely, above a critical interface re- 
sistance, the proximity-induced superconducting correlations 
in the normal metal at the Fermi level change abruptly from 
conventional even-frequency pairing to odd-frequency pair- 
ing. This result is interesting for two reasons. Firstly, the odd- 
frequency correlations penetrate much deeper into the normal 
metal region, since there is no explicit exchange field there. 
Secondly, the result provides a scenario where odd-frequency 
amplitudes are present without any interfering effects of even- 
frequency correlations. In light of the above discussion, it is 
seen that this actually resolves the two main difficulties asso- 
ciated with the experimental detection of odd-frequency cor- 
relations. 

In this work, we expand on the results provided in Ref. 39 
and address in particular three complementary issues: (i) how 
is the even- to odd-frequency conversion influenced by pair 
breaking effects near the interface, ( ii) how is the inverse prox- 
imity effect in the superconducting region influenced by the 
presence of spin-DIPS, and (Hi) how does a Fermi surface 
mismatch influence the effect under consideration? These 
questions are important from an experimental perspective, 
where non-idealities such as pair breaking effects are generi- 
cally present, and demand a numerical and self-consistent ap- 



proach. The system under consideration is shown in Fig. [T] 
The superconductor is assumed to be a conventional supercon- 
ductor such as Al or Nb, thus featuring a spin-singlet, even- 
parity (isotropic), even-frequency symmetry for the order pa- 
rameter. The insulating interface region separating the normal 
metal and the superconductor is assumed to be magnetic, e.g. 
EuO. The density of states (DOS) can be probed experimen- 
tally in various ways, for instance spectroscopically by using 
a local scanning tunneling microscopy (STM)-tip. 

This work is organized as follows. In Sec. |II A| and |III A| 
we establish the theoretical framework to be used for obtain- 
ing our results. In Sec. II B and III B we present our main re- 
sults, demonstrating that the even-odd frequency conversion 
is a robust effect which survives both in the clean and dirty 
limit, and moreover is resilient towards pair-breaking effects 
near the interface. We summarize our findings in Sec. IV 
We shall use units such that h=c=kB=l. Moreover, we use • 
for 2x2 spin-matrices, • for 4x4 matrices in Nambu-Gor'kov 
particle-hole space, and boldface notation for vectors. 

We use the quasiclassical theory of superconductivity,E^El 
where information about the physical properties of the system 
is embedded in the Green's function. For an equilibrium sit- 
uation, it suffices to consider the retarded part of the Green's 
function, here denoted g. We begin our discussion with the 
diffusive limit, after which we proceed to the ballistic case. 



II. DIFFUSIVE LIMIT 
A. Theory 

Due to the symmetry properties of g, one may parameterize 
it conveniently in the superconducting (S) and normal (N) re- 
gion in the diffusive limit.ESl Consider for concreteness an N|S 
bilayer, where we may write 
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with c = cosh(6'), s = sinh(6'), and 9 
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with Co- — cosh(f?o-). So- — sinh(6'cr). The diffusive prop- 
agators are normalized according Vq g\ — — \ where 
1 — diag(l, 1,1,1). Through this parameterization, we have 
taken into account the possibility of odd-frequency triplet cor- 
relations in the normal region, while we have employed the 
bulk solution in the superconductor. This approximation is 
valid under the assumption that the superconducting layer is 
much thicker and less disordered than the normal region, thus 
acting as a reservoir.' The gap suppression near the inter- 
face may furthermore be neglected in the tunneling limit.E^ 
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In general, the superconducting region is also influenced by 
the proximity effect, in which case a similar parameterization 
as Eq. Q is employed also in that region. We will return to 
this issue below. 

In the present case, the Green's function in the normal 
region obeys the Usadel equation 



DV(gNV.gN) + i[£P3,5n] = 0, 



(3) 



with Pa = diag(l, 1,-1,-1), and is subject to boundary con- 
ditions at the S|N {x = 0) and N|I {x = ^n) interfaces as 
follows™ 



= [.9s, 5n] 



(4) 



with f3 = diag(l, — 1, 1, — 1), at a; = and d^Oa = at 
X — d^. Here, 7 = Rb/Rn where i?B (^n) is the resistance 
of the barrier (normal region), and cJn is the width of the nor- 
mal region, while Gt is the barrier conductance. For later use, 
we define the superconducting coherence length — \/D/ A 
and Thouless energy ejh = D/d^, where D is the diffusion 
constant. Eq. Q contains an additional term compared to 
the usual non-magnetic boundary conditions in Refs. i25iand 
|26, The physical interpretation of this term is that it gives 
rise to spin-dependent phase shifts of quasiparticles being re- 
flected at the interface. Note that may be non-zero even 
if the transmission Gt — > 0, corresponding to a ferromag- 
netic insulator.l23 Later in this work, we shall also consider a 
fully self-consistent calculation where the bulk solution is not 
assumed in the superconducting region. 

Using a simplified scattering model near the interface, it is 
possible to obtain microscopic expressions for Gt and G^. 
They are related to the transmission and reflection amplitudes 
{i^*'^*, r^*'^*} on the S (N) side of the interface. For simplicity, 
we assume that the interface is characterized by N identical 
scattering channels. Under the assumption of tunnel contacts, 
one obtains from a model with a Dirac-like barrier potential 

Gt = NGqT, G^ = 2NGq (p^ - 4rVr) (5) 

upon defining T^J^a Gq = e^/{2Trh), and 

p^^lm{r^{rlr}, = Im{4(t^)*}. (6) 
The scattering coefficients take the form: 

= (fcN - - ik^Z,)/V,, tl = 2V¥k^/V^, (7) 



with file definitions ^ + + \k^Z„ , k^ = y/2rn^, 
k^ = \/2mN/iN. Here, Z„ ~ Zq + aZ^ is the spin-dependent 
barrier potential, and we define a = Z^/Zq as the polarization 
for the barrier. The ratio [G^/GtI is evaluated in Fig. |2]as a 
function of the barrier strength Zq for several values of a. We 
have used /is = /in = 5 eV and set tos = ttin to the bare 
electron mass. A Fermi-vector mismatch /is 7^ /^n between 
the materials is accounted for by an increase in Z^. As seen, 
the ratio \G^/Gt \ can be of order unity for low barrier trans- 
parencies Zq ^ 1 even for relatively weak polarizations with 
a = 10%. 
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FIG. 2: (Color online) Plot of the ratio | / Gt \ as a function of the 
barrier strength Zq for several values of the interface polarization a. 



B. Results 

We begin our analysis by discussing the weak proximity 
regime, where an analytical treatment is possible for all quasi- 
particle energies e. Thereafter, we present a self-consistent 
numerical calculation for an arbitrary proximity effect, incor- 
porating pair-breaking mechanisms and the depletion of the 
superconducting order parameter near the interface region. In 
the linearized treatment, one assumes that the deviation from 
the bulk Green's function in the ferromagnetic region is small. 
This permits us to write the retarded Green's function on the 
form 



Here, we have defined 
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Under the assumption of an equilibrium situation, the Keldysh 
Green's function is given by 



=[r-.r]tanh(/3£/2), 



(10) 



where /3 ~ 1/T is inverse temperature. The advanced compo- 
nent is = —{p^g^p^y . The linearized Usadel equation^ 
may be written as 



Ddlf± + 2ie/± = 0, 



(11) 



and is to be supplemented with the boundary condition ob- 
tained from Eq. Q 
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at X = while dxf± = at x = d^- Here, f± = /t ± 
/s where ft is the 5^ = triplet component and fs is the 
singlet component of the anomalous Green's function. Since 
the diffusive limit is considered, the singlet component has an 
even-frequency symmetry while the triplet component has an 
odd-frequency symmetry. The odd-frequency component has 
previously been predicted to appear in S|F layers, but we now 
show that the presence of a magnetically active barrier region 
induces an odd-frequency component in a S|N layer, with a 
much longer penetration depth. We find that the solution for 
the Green's function reads 



\G<b\ > Gi 



f± 



-^gUik{x-2dti) _|_ g-ifea:l 



ik-fdt^{l - c-2ife*) + (c ± iG^/Gt){1 + e-2iferfN) ■ 



(13) 



Here, k = y/2iejD. For a spin-inactive barrier, G^ = 0, 
we obtain /+ = — /-, such that ft — 0. However, the pres- 
ence of Gcf, induces the odd-frequency component in the nor- 
mal layer. The decay length here is not dictated by the mag- 
netic coherence length ~ \J~Djh as in an S|F layer, but 
by — \/D/e as in an S|N layer. This allows the odd- 
frequency component to penetrate much deeper into the N 
layer than into the F layer. The simplest experimental man- 
ifestation of the odd-frequency compone nt is pro bably a zero- 
energy peak in the local density of states .GMEllnSlF layers, 
where this phenomenon has been discussed previously, a clear 
zero-energy peak is unfortunately often masked by the simul- 
taneous presence of singlet correlations {fs), which tend to 
suppress the density of states at low energies. In the present 
case of a spin-active interface in an S|N junction, however, 
Eq. (13 1 suggests a remarkable effect. Consider e = 0, for 



which fc = 0, s = 1, and c = 0, leading to the result 

f± = Gt/ G^ 



(14) 



under the assumption that ^ 0. This equation conveys 
a powerful message, namely that at the Fermi level, the sin- 
glet component is absent while the triplet component remains. 
Moreover, the latter is determined simply by the ratio of Gt 
and G^. Consequently, this should provide ideal circum- 
stances for direct observation of the odd-frequency compo- 
nent, manifested as a zero-energy peak in the local density of 
states. 

So far, we have limited ourselves to the weak proximity 
effect regime. We now consider an arbitrarily large proximity 
effect. In this case, the Usadel equation reads 



Ddle„ + 2iesmhe^ = 0, 



(15) 



while the boundary conditions become 



Gd 



^dt^dxOa ^ (cSa ~ CFSCa) + ICr-p^Sc (16) 

Lit 

at a; = and dxda = at a; = cJn- A general analytical 
solution of the above equation can hardly be obtained, but it 
may be solved at zero energy. For e = we find pairing 
amplitudes that are either purely (odd-frequency) triplet for 



/,(0)=0, /t(0) 



Gt ■ sgn{G^) 



Gi 



(17) 



or purely (even-frequency) singlet for jG^I < Gt, 
1-Gt 



/s(0) = 



/t(0)-0. 



(18) 



Thus, the presence of G^ induces an odd-frequency compo- 
nent in the normal layer. The remarkable aspect of Eqs. (17 1 



and (1 8 1 is that they are valid for any value of the width 
below the inelastic scattering length, and for any interface pa- 
rameter 7. Thus, the vanishing of the singlet component is a 
robust feature in S |N structures with spin-active interfaces, as 
long as jG^I > Gt- Without loss of generality, we focus on 
positive values of G^ from now on. The DOS is given as 



N{e)/No = J2Mc.}/2, 



(19) 



thus yielding 



N{e = 0) 



Re 



Gd 



y/G} - Gj 



(20) 



At zero-energy, the DOS vanishes when G^ < Gt, which 
means that the usual minigap in S|N structures survives. How- 
ever, the zero-energy DOS is enhanced for G^ > Gt since the 
singlet component vanishes there. 

We suggest the following qualitative explanation for the 
mechanism behind the conversion between even- and odd- 
frequency correlations. The superconductor induces a mini- 
gap oc Gt in the normal metal, while the spin-active barrier 
induces an effective exchange field oc G^. The situation in 
the normal metal then resembles that of a thin-film conven- 
tional superconductor in the presence of an in-plane external 
magnetic field,^^ with the role of the gap and field played by 
Gt and G^, respectively. In that case, it is known that super- 
conductivity is destroyed above the Clogston-Chandrasekhar 
limit,^^ as the spin-singlet Cooper-pairs break up. In the prox- 
imity structure we consider here. Cooper-pairs persist above 
this limit as they are induced from the superconducting region 
where the exchange field is absent. However, these Cooper 
pairs are modified strongly by multiple scattering from the 
spin-active interface, and above a critical ratio G^ /Gt — 1 
spin-singlet pairing is no longer possible in the N region at 
the chemical potential. It is then replaced by spin-triplet pair- 
ing, which must be odd in frequency due to the isotropization 
of the correlation in the diffusive limit. We observe coexis- 
tence of the exchange field and spin-singlet even-frequency 
superconductivity as long as G^ is below the critical value of 
G^ = Gt- At the critical point, the DOS varies as 1/ -\/[ef and 
diverges at e = 0. Thus, we find that there is a natural sep- 
aration between even-frequency and odd-frequency pairing in 
the normal metal at a critical value of the effective exchange 
field G0. This agrees with the interpretation of G^ in Ref.|29l 
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as an effective proximity-induced exchange field. Note that 
the above expressions are valid also for 0: we obtain 

fs = 1 and /t = as demanded by consistency. 

The full energy-dependence of the DOS may only be ob- 
tained numerically. In addition, it is of interest to see how 
robust the predicted even- to odd-frequency conversion is to- 
wards the inevitable depletion of the superconducting order 
parameter near the interface in addition to non-ideal effects 
such as the presence of inelastic scattering. To investigate 
this, we solve the Usadel equation and the gap equation self- 
consistently in both the normal and superconducting region. 
Since we are no longer considering the bulk solution in the 
superconducting region, it becomes necessary to specify the 
width ds of the superconducting layer, the spin-dependent 
phase shifts on the superconducting side of the interface, 
and also the bulk resistance i?s of the superconductor. The 
Usadel equation on the N side satisfies Eq. (j3]l, whereas on 
the S side an additional term A is added inside the commuta- 
tor in the second term of Eq. (j3]l. Inclusion of spin-orbit cou- 
pling effects may be done similarly by including a term CTso 
(see Ref. 45 for a detailed treatment and expressions for such 
terms). The superconducting order parameter is determined 
self-consistently by solving the Usadel equation in conjunc- 
tion with the gap equation: 



A = 



iVpA 



/ d£tanh(;9£/2) VcrRe{sinh(6'^)}, (21) 

JO 



where we choose the weak coupling-constant and cut-off en- 
ergy to be N-pX = 0.2 and oj/Aq = 75. Within our numerical 
scheme, self-consistency is typically achieved after 10 itera- 
tions. We account for inelastic scattering by the parameter 
d/Ao = 10^3, where e + i5. 

The diffusion coefficients Z^n and Ds are in general dif- 
ferent. At the S|N interface, the boundary condition on the 
normal side now reads: 



Rb 



while on the superconducting side, one has: 



its (-fT 



G5 



(22) 



(23) 



The magnitude G^ of the phase-shifts induced in the super- 
conducting region are equal to G^ in the absence of a Fermi- 
vector mismatch, but will in general be different. The normal- 
state conductivities are given by 



"N(S) 
-Rn(S)^ 



(24) 



with A as the interface area and i?N(S) is the normal-state resis- 
tance. Since it is reasonable to assume that the barrier region 
features a higher electrical resistance than the bulk of the ma- 
terials, we shall set Rq/R^ = 4 in what follows. Moreover, 
we fix the width of the normal layer to c^n/^s = 1-0. 

Due to an inverse proximity effect, the superconductor 
should also be influenced by the presence of G^ 7^ 0, and one 



expects that an odd-frequency triplet component would be in- 
duced near the interface on the superconducting side. There- 
fore, we will also study how this inverse proximity effect is 
manifested in the superconducting DOS. We will focus on the 
influence of the spin-DIPS G^, considering an equal magni- 
tude of spin-DIPS in both regions, i.e. G^ = G^. Consider 
first a situation where the superconducting region acts as a 
reservoir and is very weakly affected by the proximity effect. 
To this end, we set ds/Cs — 5.0 and ctn/cs = 0.2, ensuring 
in this way that both ds 3> and that the superconducting 
region is less disordered than the normal region. 

The results are shown in the top row of Fig. |3] where we 
plot the DOS in the superconducting region, the normal metal 
region, and also the spatial depletion of the order parameter 
The DOS is plotted at x — —da, in the superconducting re- 
gion and X — dj^ in the normal metal region, and may be 
probed by tunneUng spectroscopy measurements through an 
insulator In the superconducting region, the results are vir- 
tually independent of G^ in the present case of a reservoir 
modeled by ds/^s = 5.0, so we consider only G^ = there. 
As seen, both the inverse proximity effect and the gap deple- 
tion are negligible. However, the DOS in the normal metal 
region is highly sensitive to the presence of G^. In particular, 
the low-energy DOS displays a strong dependence on the ratio 
G^/Gt- We will comment further on this below. 

In the bottom row of Fig. [3] we investigate a scenario where 
the superconducting region no longer acts as a reservoir, and 
where the proximity effect is expected to be substantial in both 
regions. To this end, we fix ds/^s = 1-0 and o'n/o's = 1-0. In 
this case, the proximity effect in the superconducting region is 
much stronger than in the reservoir case of ds/^s — 5.0, and 
the depletion of the superconducting order parameter is more 
pronounced. In particular, the DOS at Fermi level is no longer 
zero and depends on the value of G^. However, both the DOS 
and the superconducting order parameter remain quite insen- 
sitive to a variation in G^. In the normal metal region, the 
behavior is similar to the reservoir case, although the peak 
structure at zero-energy now appears for a lower value of G^. 

We are particularly interested in seeing if the even- to odd- 
frequency conversion predicted from the analytical treatment 
in Sec. |Il] is equally pronounced in this numerical, self- 
consistent treatment. To this end, we plot in Fig. |4]the zero- 
energy DOS in the normal metal at x = dN as a function 
of Gcf, I Gt for both the case of a superconducting reservoir 
{ds/£.s = 5.0, ctn/cts = 0.2) and a thin layer (ds/Cs = 
1.0, (Tn/cs = 1-0). In both cases, the transition from a fully 
suppressed low-energy DOS to an enhanced low-energy DOS 
appears at 



Gcf, — tiGt, 



(25) 



where i] < 1. This is a clear signature of the transition 
from pure even- to pure odd-frequency correlations. The cor- 
responding behavior of the anomalous Green's function is 
shown in Fig. |5] where we have included inelastic scatter- 
ing and solved self-consistently for the superconducting orer 
parameter As seen, the correlations undergo a rapid transition 
from singlet to triplet at G^/Gt = tj, with ?/ e {0, 1}. 

From our above findings, it then follows that the even- to 
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Superconductor at x = —d^ Normal metal at x = dj^ 
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FIG. 3: (Color online) (a) and (c): Plot of the DOS in the superconductor at a; = —ds- (b) and (d): the DOS in the normal metal at a; = d^. 
Insets: the spatial profile of the superconducting order parameter. The black arrows indicate an increasing value of / Gt- In the top row, we 
model a scenario where the superconductor acts as a reservoir, we have set ds/Cs = 5.0, dN/6 = 1-0, a^^/as = 0.2, and G| = = G<^. 
In the bottom row, we model a scenario where the proximity effect is expected to be substantial in both the N and S regions, we have set 
ds/Cs = 1.0, dN/Cs = 1.0, rJN/as = 1.0, and G^ = G^ = G4,. 
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odd-frequency conversion persists also for thin superconduct- 
ing layers comparable in size to the coherence length ^s, as 
long as the inverse proximity effect is relatively weak, with a 
concomitant shift in the critical interface resistance where the 
pairing transition occurs. 



III. BALLISTIC LIMIT 



A. Theory 



FIG. 4: (Color online) Plot of the zero-energy DOS in the nor- 
mal metal versus G^/Gt for the two cases of ds/^s = 5.0 (with 
ctn/cts = 0.2) and ds/^s = 1-0 (with at^/as, — 1.0). As seen, an 
abrupt transition occurs at a value G,f, = tiGt, where r? < 1. 



Turning our attention now to the ballistic limit, our strat- 
egy will be to solve the Eilenberger equation and supple- 
ment the solution with boundary conditions obtained by 
means of th e tS-m atrix method elaborated upon in a num- 
ber of works .'^ES^^'l^ The retarded Green's function g = 
is in this case most conveniently parameterized by Riccati- 
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FIG. 5: (Color online) Plot of the anomalous Green's function at 
zero energy (Fermi level) as a function of G^/Gt at a; = dN. In 
(a) we use ds/Cs = 5.0, (Tn/o"s ~ 0.2, while in (b) we use ds/Cs = 
1.0,(7n/(Ts = 1.0. As seen, in both cases a transition occurs from 
singlet to triplet coiTelations at G^i/Gt = rj where rj £ {0, 1}. We 
have included inelastic scattering and solved self-consistently for the 
order parameter. 



amplitude J^^I^^Ell {^^^ x^^^ where 

'A/; (1+77) 2A/;7 
-2Uj_ (1+77) 

and the normalization matrices read: 

7y: = ( l-77)-i, = ( 1-77)"^ 



(26) 



(27) 



Here, we use the notation of Ref. 24 assuming the Green's 
function to be normalized as 5^ — — tt^I. The Eilenberger 
equation for the propagator in the normal region, g = g^, 
reads 



(28) 



where = diag(l, 1,-1,-1). For the boundary conditions 
at the interface we closely follow the 5-matrix approach in 
the form presented in Ref. [38| The scattering approach de- 
scribes the system by separating it into a scattering region, 
which cannot be described within quasiclassical (QC) theory, 
and asymptotic regions on both sides of the interface, where 
QC theory is applicable. The scattering region must be small 
compared to the coherence length. It must also extend far 
enough into the asymptotic region, such that the QC theory is 
applicable. The 5-matrix approach essentially consists of de- 
termining the unknown Riccati amplitudes corresponding to 
trajectories starting at the interface and moving into the bulk 
on each side by relating them to the known Riccati amplitudes 
describing trajectories starting in the bulk and moving towards 
the interface. These two sets of amplitudes are related pre- 
cisely via the iS-matrix. 

The details of the iS-matrix depend on what kind of in- 
terface is considered. For our purposes, we shall consider a 
quite general model. Namely, an interface which is (i) par- 
tially transmitting (non-ideal), ( //) specular (parallel momen- 
tum is conserved), and (Hi) spin-active (giving rise to spin- 
mixing and spin-filtering effects). The 5-matrix is evaluated 



at the Fermi level in the quasiclassical approximation, i.e. 
S = S{pp), and can be written as 



S = 



TnS — i?N 



(29) 



The indices S and N refer to the superconducting and nor- 
mal metallic side of the interface, respectively. Thus, Sss 
describes reflection processes at the superconducting side of 
the interface, whereas Ssu describes transmission from the 
superconductor to the normal metal. The elements Sij with 
{i, j} E {S, N} are diagonal in particle-hole space according 
to 











where py denotes the component of the momentum parallel to 
the interface, and the superscript tr denotes matrix transpose. 
In the presence of an inversion symmetry within the interface 
plane, the sign of p|| is unimportant. In general, interface scat- 
tering may allow for spin-flip processes, namely when spin- 
rotation invariance is completely broken in the system under 
consideration. The details will depend on the micromagnetic 
properties of the interface. Here, we will treat the common 
case that spin-rotation invariance is only partially broken, i.e. 
it is still present with respect to rotations around the axis along 
the magnetic moment of the interface. Choosing our quanti- 
zation axis along this direction, the scattering matrix is also 
diagonal in spin-space and has the general form: 











(30) 



Current conservation requires unitarity of the scattering ma- 
trix, i.e. the parameters defined by this equation are not inde- 
pendent. Moreover, the physical results obtained from qua- 
siclassical theory must be gauge invariant in the following 
sense. We may transform the 5-matrix by: 



5' = 



e"'2/2i_ 



e"'2/2l 



(31) 



without changing the solutions of the quasiclassical boundary 
conditions. This additional gauge freedom is related to the 
fact that only the envelope of the wave function enters quasi- 
classical quantities sufficiently far away from the interface. A 
transformation according to Eq. ( (3T| only changes the wave- 
function on either side of the interface by a scalar phase factor 
and thus is irrelevant on the quasiclassical level. The same 
gauge freedom can be used to show that the precise definition 
of which part of the system is to be included in the scatter- 
ing region (within the abovementioned restrictions) does not 
influence any physical quantity calculated within QC theory 
(using a general form of the 5-matrix, for the current problem 
this is shown in Ref. 38). Exploiting unitarity and the above 
gauge freedom, we arrive at the following parameterization of 
the 5-matrix 



S = 



r e 



ii}sgij2 ^ gi(i5sN 0:^+0' l)/2 



^ gi(i5NS o:,-<^' l)/2 



-r e 



(32) 
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a; = 



X — djv 



FIG. 6: (Color online) Illustration of the incoming and outgoing 
Riccati-amplitudes in the NjS bilayer. Lower-case amplitudes (7, 7) 
should be integrated towards the interface at a; = 0, whereas upper- 
case amplitudes (F, T) should be integrated away from the interface 
at X = 0. In the superconductor, we use the bulk solution. All am- 
plitudes are homogeneous in direction parallel to the interface. 



with r = diag[r.|-, r^] and t = diag[t.|-, t^], and is the third 
spin PauH -matrix. Above, (p' arises due to a possible contribu- 
tion from a vector potential when a magnetic field is present in 
the interface region. This contribution is independent of spin, 
and originates from the time-reversal symmetry breaking by 
the magnetic field at the interface. It gives an extra phase to 
the anomalous components, and basically corresponds to the 
magnetic flux through the interface cross section. It is irrel- 
evant for our purposes, i.e. the behavior of the DOS, but we 
have kept it for the sake of generaUty. Unitarity requires 



tl + rl 



i9 



NS ■ 



t9. 



SN 



1, 



(33) 
(34) 



which implies 6 free parameters. 

With the 5-matrix in hand, the remaining step is to write 
down the appropriate boundary conditions which serve as the 
link between the incoming Riccati amplitudes 

In = In^^'II'"^^''^'^)' In = 7N(P|hPx,e,a;), (35) 
and the outgoing Riccati amplitudes 

En = I^^{P\\,Px,e,x), f^ = V^{p\\,-p^,e,x). (36) 

measured with respect to the S|N interface (for the notation 
see Fig. [6]). The general solution of the Eilenberger equation 
in the normal metal region reads 



= 7jdN)e-2-(-*)/"- 



(37) 
(38) 



for the trajectories along 6, whereas for trajectories along vr— 6* 
we obtain 

lN(^)=2N(^N)e-2'^(^-'^^)/''-, (39) 

rN(a;) = ^N(o)e2'^"/''^^ (40) 

where we defined v^x ~ vy^cosO, and — 7r/2 < 6 < n /2 
is assumed. Here, and in the following, we suppress the pa- 
rameters px and e in the argument list. All amplitudes are 



independent of p| | . The bulk solution is used for the incoming 
Riccati amplitudes on the SC-side, 



Is 



-Is 



Ao 



(41) 



where we used a real gauge for the superconducting order pa- 
rameter Aq, and is the second spin Pauli-matrix. As shown 
in Ref.jl^ the foUowing boundary conditions at a; = hold: 



■SN 

1 



EN(0)=7NN + rN^s7s(0)7U (42) 

rN^s-i;,s [i^IsWl^s 

i;.=E^.a.(o)5,. 

where j, k and I run over {N, S}. Analogous equations hold 
for rN(0). At a; = cJn we assume perfect and non-spin-active 
reflection, hence the boundary conditions are trivial: 

1^{dn) = I.^{dn) = tnidn) = %{dN) = (43) 

and result in the following relations between amplitudes at 

X = and x = d^: 



EnIO) 



736 



7,(0) 



2ied^/vp 



(44) 
(45) 



Replacing En(0) ™d 7n^^-' ™ 63 according to this re- 
lation yields a quadratic equation in 78,0- whose solutions can 
be determined analytically.'^ 



B. Results 

The odd-even frequency conversion which was shown to 
take place in the diffusive limit also occurs in the ballistic 
limit, as we show in the following. In this case, we obtain 
the retarded Green's function using the formalism described 



in Sec. Ill A There we derived an equation from Eqs. ( |42| ) 
and ( |44] i for the Riccati amplitudes in the normal metal region 
that determine the proximity amplitudes. Following Ref.[38l 
we obtain analytical expressions for the anomalous Green's 
function in the N region. The energy-resolved DOS at the 
outer boundary of the normal layer can then directly be ob- 
tained from the Riccati amplitudes via: 

- -Im-(.gB> = y ((I-IbIb) '(I + IbIb)) 



No 



where (•) denotes the Fermi-surface average given by: 



1 



d^p'p 



No J PS (27rfi)3|^F(p^)|^'^' 
with the local density of states in the normal state, 

d^p'p 



No 



PS {2nhr\vFip'^)\- 



(46) 



(47) 



(48) 
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It is important to realize that only the singlet and the 5*2 = 
triplet component (in a basis where the z-axis is along the 
quantization axis) will be induced in the normal part of the 
system, since the magnetization of the barrier is uniaxial and 
has no inhomogeneity. We may then write 



7p 



7+ 

-7_ 



7+ 
7_ 



(49) 



which can be inserted into Eq. (42 1, and similarly for 7, 



Using the scattering matrix definea in Eq. ( |32| i and focus 
ing on subgap energies |e| < Aq, where 7s = le'* with 
^I^ = arcsin(e/Ao) (here -7r/2 < * < ixjl), we get two 
decoupled equations for ^„ : 

2u, , 



1 = 0, 



(50) 



where u G {- 

Ua-{e) = sin 



}, and the function Ua{e) is defined as 



+ r^ri sin 



/2£dN 
V VFx 



(51) 

and the variable 



Here, we have defined 'd± = ^{'^n ± 
a is to be understood as a factor ±1 for cr = ±. Eq. (50 1 is 
solved by: 



7^ 



± 



1 



i) 



(52) 



We can write down an equation analogous to Eq. (50i for 70-, 
2m„ 



2i4>' 



1 = 0, 



with u„{£) = Uc{—e). Noting that u-{—e) = —u+{e), it 
follows that = — 7+G^"^ . The correct sign is obtained by 
requiring (i) that the symmetry relation holds between 7^ and 
7cr, and (ii) that the momentum and spin resolved density of 
states, e.g. 



No 



1 



■ 7+7- 



2 ^2ld 



1 + 7+7- 



1 + 7+e- 
1 - 72 e2"?; 



(54) 



must be positive. Also, we must demand 7cr — > when 
t^ti — > 0, as in that case the two regions become completely 
decoupled and the proximity effect should be zero. It follows, 
that the appropriate solution is 



ja^e sgn[Mo 




(55) 



Having obtained the correct solution for the 7^ and ja- 
quantities, we turn to the anomalous Green's functions. From 



the parameterization of the Green's function Eq. (26i, we 
identify 

fa = -27ri ^\ = -27ri , (56) 



where we have defined f± = fs ± ft- Inserting the solutions 
of 7ct and 7^ from above, we obtain 



fa{e) = ITT 



t^t^c-'^'sgn[u,(£)] 



(57) 



At this point, it is instructive to consider the angular de- 
pendence of the Green's function. From the boundary condi- 
tions in Sec. |III A| and the general form of 7^^ and Fj^, it is 
^^((iN) for all values of 6. In effect, this 
although we used 



-b; seen that 7 (di 



means that Eq. ( 57 1 is valid for any value 



e.g. 7j^(dN) to obtain them, which only is defined in the range 
— 7r/2 < 6* < 7r/2. It is also worth to note that at the outer sur- 
face of the bilayer, all singlet components are even-frequency 
while all triplet components are odd-frequency. While this can 
be shown analytically from the above equations, one may also 
understand it intuitively from the fact that one has specular re- 
flection at the outer surface such that all components must be 
even in momentum there. 

We obtain the energy-resolved DOS at the outer boundary 
of the normal layer from the equations presented above. 



N{e) 

Nn 



(58) 



and the pairing amplitudes from Eq. ( [ST] ). In order to inves- 
tigate the even- to odd-frequency conversion at the chemical 
potential, we are in particular interested in their value at e = 0. 
We obtain for |ito| > i^^^- 

/3(s = 0)=0, ftie^O)^..^^^^, (59) 



(53) 

whereas for |uo| < t'fti- 



fs{£ = 0) = TT 



/t(e = 0)=0. (60) 



Here, the parameter uq is given by 

Mo = sni + r^r; sm ) . (61) 



In the case that both and i9s are of order of t^t^, and the 
system is at the same time in the tunneling limit, we can ex- 
pand all quantities up to (t^t^)^ and thus recover the results 
of Ref. 39 that uq in the above expressions for the pair ampli- 
tudes and for the density of states can be replaced by ?9n- Note 
that the scalar phase 0' was set to zero in Ref. |39l as it has no 
consequence for the behavior of the DOS. 

Considering a realistic interface, it is clear that both the 
transmission coefficients and the spin-mixing angle i^n de- 
pend on the angle of incidence 9. A systematic study of 
this angular dependence was performed in Ref 21 While 
the expressions Eqs. ( |59l ) and ( |60l ) are valid even for angle- 
dependent quantities and dj^, the Fermi-surface average 
can in fact add the possibility of a simultaneous presence of 
both triplet and singlet correlations at zero energy originat- 
ing from different incidence angles 9; i.e. some trajectories 
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FIG. 7: (color online) Definition of the parameters for the spin- 
dependent barrier. For Figs. |8|lT| we use V5 = 0.2i?F, and di = 

2/fcF,S. 



may contribute to the singlet component while others con- 
tribute to the triplet. In order to discuss the conditions under 
which that happens, we performed calculations for an inter- 
face layer modeled by a spin-active box-shaped potential of 
width di, with di ^ Ap. Its height is spin-dependent and given 
by Ui\ = V]_, Uii — Vi + Ji, where Ji is the interface exchange 
field (see Fig.|7]for the notation). 

For all results presented below, the lower of the two poten- 
tial barriers is Vi = 0.2i?F, and the interface width is dy = 
2AF,s/27r. We assume for simpUcity equal band masses all 
over the system, and isotropic Fermi surfaces. Thus, the en- 
ergy dispersions are in the superconductor k^/2m, in the nor- 
mal metal E^+k^ /2m, and in the barrier E^+Ui^^^^+k^ /2m, 
where i?p = fcp s/2'7^■ The constant determines the Fermi 
surface mismatch between the superconductor and the normal 
metal, with Fermi wavevectors fcp.s and fcp.N for the supercon- 
ductor and the normal metal, respectively. For such a model, 
the parameters in Eq. (|29]l are given by 



-SN 



1 - 



1 - 



NS' 



(62) 
(63) 



where p. with j e {S,N}, and u are diagonal spin-matrices 

Withpj^cro- = (fcj -lKcr)/(fcj +lKcr), andWcrcr 



Here, k„ = [fc| gi/i.^/E'F + k'^\Y''^, and % - ^apj 



: exp(-KCTC?i). 

We now turn to the discussion of our results, shown in 
Figs.[8pT| The calculations were obtained for various Fermi 
surface geometries. In Fig. |8]we present results for the case 
^F,s ^ ^F,N, i-S- when the Fermi surface mismatch is large, 
and the Fermi surface in the superconductor is much larger 
than that in the normal metal [see inset in Fig. [8|a)]. A 
more moderate mismatch is assumed in Fig. |9] with fcp.N — 
O.Sfcp.s. In Fig. 10 we consider the special case of no Fermi 
surface mismatch, i.e. fcp s = fcp.N- Finally, in Fig. 11 
we consider the case opposite to Fig. [8] namely a strong 
Fermi surface mismatch where the Fermi surface in the nor- 
mal metal is much larger than that in the superconductor, 
^F,s ^ ^F.N- In each figure, we present in (a) the energy 
resolved DOS for several spin polarizations of the interface 
barrier. The energy scale of interest is the Thouless energy 



of the normal metal layer, exh = fivp,]<i/2d'^, times the pair 
transmission amplitude from the superconductor to the nor- 
mal metal, t^t^. This quantity depends on the impact angle 
9; for definiteness we use as energy scale Eq ~ t^ti^Th for 
normal impact. In (b) we show the quantity t^^t^ —\u[)\, where 
Mo = u+(£ = 0) = — u_(£ = 0) is defined in Eq. ( |6T| l. The 
plotted quantity controls the transition from even-frequency 
singlet to odd-frequency triplet correlations at the chemical 
potential (e = 0), as will be shown below. The parameters 
z9s,N and r^^^ depend on fc||, the momentum component paral- 
lel to the interface, which is conserved in the scattering pro- 
cess. Consequently, the parameter t^t^ — |uo| depends on fc|| 
as well, and we show in the figure this dependence. In (c) we 
show the even-frequency singlet and the odd-frequency triplet 
superconducting amplitudes at the chemical potential [we plot 
the real quantities /s(e) = and i/t(e — 0))] at the outer sur- 
face of the normal metal. Finally, in (d) we show the local 
DOS at the chemical potential normalized to the normal state 
DOS, N{e — 0)/No, again at the outer surface of the normal 
metal. 

We proceed with the discussion of the results. We recall 
first the known behavior for zero interface spin polarization, 
Ji = in Figs. 8pT When all trajectories in the normal 
metal are partially transmissive (Figs. 8][T0i, the DOS is zero 
at the chemical potential, e = 0, and shows an increase to fi- 
nite values as function of energy e. This increase is directly 
associated with the behavior of the topmost ( Ji = 0) curves 
in Figs. 8pO b) for glancing impact, fc||/fcp,N ~ 1- When 
there are non-transmissive trajectories present in the normal 
layer (Fig. Ill, the DOS at the chemical potential is finite. 
In this case, as fcp s < fcp,N, there is a background DOS re- 
sulting from the non-transmissive directions, fc|| > fcp,s, in 
the normal metal; this background contribution is not associ- 
ated with any superconducting pair correlations, and is nearly 
constant in energy and nearly temperature independent (con- 
sidering typical superconducting energy scales). All changes 
of the DOS related to superconductivity take place on top of 
that background contribution. 

When the spin polarization of the interface increases to non- 
zero values, we can define three characteristically different re- 
gions of interface spin polarization Jj. We turn our attention 
to panels (b) of Figs.[8][TT] which show the quantity t^t^— |uo| 
as a function of fc||. For directions where this quantity is pos- 
itive, according to Eq. (|59]l pure singlet correlations are cre- 



ated at the chemical potential in the normal metal, whereas 
for directions where this quantity is negative, according to 
Eq. ( |60l ) pure odd-frequency triplet correlations are created 
at the chemical potential in the normal metal. We can clas- 
sify the curves into three groups, depending on the value of 
Ji. We first have a region where t^t^ — \uo\ is positive for 
all fc|| (region I; e.g. Ji = 0.2£'p in Fig.lSll; second a region 
where t^t^ — \uo\ is positive for some, and negative for other 
values of fc|| (region II; e.g. Jj = 0.2_Ep in Fig. 10 1; and third 
a region where t^t^ — \uo\ is negative for all fc|| (region III; 
e.g. Ji > 0.3Ep in all four figures). 



In Figs. 8|TT c) we show the singlet (/g) and triplet (ft) 
component of the momentum-averaged, i.e. s-wave, corre- 
lation functions at the chemical potential. In general, there 
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FIG. 8: (Color online) (a) Energy-resolved DOS in the normal metal for different values of the interface exchange field Jj. The energy scale is 
eo = [it^i^Th] 1=0), with the Thouless energy eTh ~ fiwF,N/2dN- (b) Interface parameter tfti — \uo \ as a function of trajectory impact 
(parameterized by fcj|). (c) Singlet and triplet component of the anomalous Green's function at e = as a function of Ji. (d) Density of states 
at the Fermi level, A''(0). The Fermi surface mismatch is fcp,N ~ O.lfcp^s- The inset in the lower left corner of panel (a) is meant to illustrate 
the Fermi-surface mismatch. In (a)-(d), the interlayer thickness is di — 2Ap/27r, and the interface potential Vi — 0.2£'p. The width of the 
normal layer is = fiup,N/Ao. 
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FIG. 1 1 : (Color online) The same as in Fig.[8] however with a different Fermi surface mismatch, fcp 



lOfep' 



also exist higher order even-parity components which behave 
qualitatively similar When increasing Jj in region 1, it can 
be seen from Fig. [8] and |9] that the pair correlations at the 
chemical potential stay purely singlet, and the DOS at the 
Fermi level, shown in (d), stays zero. When Ji enters re- 
gion II, there is a strong mixing between singlet and triplet 
amplitudes, and the DOS at the Fermi level rises to non-zero 
values. Finally, when Ji is above J^h (region 111) the singlet 
correlations vanish identically at the chemical potential, and 
pure odd-frequency triplet amplitudes remain, when the DOS 
is larger than its normal state value. The transition from the 
region III can be identified as a sharp decrease of the DOS as 
function of Jj from a maximum value in Figs.[8pT| 

Region I only exists for /cf,n < kp,s (Figs. [8]and|9jl. As 
seen from Fig. [Sjd), as long as /cf,n ^ ^f.s the DOS shows 
just as in the diffusive case a rather sharp transition from 
A^(0) = to a value above the normal-state DOS as a func- 
tion of the interface exchange splitting Jj. The existence of 
region II in Figs.[9pT|is due to the fact that the mixing angle 
drops slower with impact angle than the transmission. It is 
characterized by a zero crossing of the parameter t^t^ — |mo| 
as function of parallel momentum fc| | . For increasing fcp.N re- 
gion II extends to lower values of Ji, and when fcp n > fcp.s^ 
region II starts at Jj = and extends to a critical value J^^a ■ 
This is due to the fact that for any small Ji ^ there are neg- 
ative values of t^ti — \uq \ for the largest transmissive fc||. For 
fcp.N > fcp.s this can be understood easily because the trans- 
mission probability goes to zero whereas the spin-mixing an- 
gles stay finite when fc|| approaches fcp.s. For Ji > JcHt the 
system is in region III. For any mismatch between the Fermi 
surfaces, there is a critical value Jcrit- 

It is interesting to note that, although both the spin-mixing 
angles and the transmission probabilities vary with Fermi sur- 
face mismatch, for a box-shaped potential the critical value 
Jcrit does not depend on the ratio /cF.N/fcp.s- Thus, it has the 
same value, Jcrit ~ O.SEp, in Figs 
mined by the condition that t^t^ 
Eqs. (|62]i and (|63| into Eq 



This value is deter- 
for fc| I =0. Inserting 



Uq 

(|29|, and using Eqs. ([32]i and (|6T), 
this condition leads to the following implicit equation for the 



value of Jc, 



4=1 



1 



1 



1 



{vi - ly^) sinh [{vi + v^)Si] 

+ i^t) sinh [in ~ VA^)^{\ , (64) 



with the parameters ~ (Vj + Jcrit)/ v\ — ^/Vi/Ep, 
and 6i ~ fcp g di. Solutions of this equation are shown in 
Fig. 12 We find that the transition occurs earlier for thicker 
interfaces. This is because the transmission decreases with in- 
terface width, while the mixing angle is actually enhanced to 
some extent, as discussed in Ref.|2T| In order to both achieve 
a satisfying transmission, and to have realistic values for the 
exchange field, difcp.s should be between 1 and 2. The re- 
markable robustness of the critical interface spin-polarization 
with respect to the Fermi-surface mismatch might simplify the 
experimental task to observe this effect, as the usual restric- 
tions for finding suitable materials to match at the interface 
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FIG. 12: (Color online) The critical value Jcrit as function of in- 
terface thickness for various strengths of the interface potential. The 
curves were obtained by finding numerically the solution of Eq. l|64^. 
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are relaxed. We caution, however, that the above strict inde- 
pendence on the Fermi-surface mismatch might be relaxed for 
more realistic interface potentials. 

IV. SUMMARY 

In this work, we have provided a comprehensive treatment 
of the proximity effect in a system consisting of a normal 
metal (e.g. Cu) in contact with a conventional s-wave super- 
conductor (e.g. Al) through a spin-active interface. Such a 
spin-active interface is incorporated by using, e.g., a ferro- 
magnetic insulator such as EuO. We have shown that based 
on the self-consistent calculation in the diffusive regime, the 
even-odd frequency conversion first predicted in Ref.j^is ro- 
bust even when taking into account pair-breaking effects near 
the interface which cause a depletion of the superconducting 
order parameter. Although the conversion relies crucially on 



interface properties which vary considerably with the impact 
angle of incident quasiparticles, it is generically robust against 
Fermi-surface averaging in the clean limit. Moreover, we 
show that the conversion takes place even when the supercon- 
ducting region does not act as a reservoir, i.e. when the thick- 
nesses of the superconducting and normal layers are compa- 
rable. Our findings suggest a robust and simple method of 
obtaining a clear-cut experimental signature of odd-frequency 
superconducting correlations. 
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